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Abstract 

We establish a correspondence between the evolution of the distribution of 
eigenvalues of a x matrix subject to a random Gaussian perturbing ma- 
trix, and a Fokker-Planck equation postulated by Dyson. Within this model, 
we prove the equivalence conjectured by Altshuler et al between the space-time 
correlations of the Sutherland-Calogero-Moser system in the thermodynamic 
limit and a set of two- variable correlations for disordered quantum systems 
calculated by them. Multiple variable correlation functions are, however, 
shown to be inequivalent for the two cases. 
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In a series of recent papers, an interesting generalization of the problem of eigenvalue 
statistics for complex quantum systems has been introduced by Altshuler, Simons and 
coworkers They consider the change in the positions of the eigenvalues in response 
to an external perturbation; after an appropriate normalization of the perturbing potential, 
they show that the evolution of the eigenvalues of the system as a function of the strength of 
the perturbation is universal. By treating this extra parameter (the strength of the perturba- 
tion) as a time-like coordinate, they conjecture that this problem is identical to the ground 
state dynamics of an integrable one- dimensional interacting many body quantum model, 
the Sutherland model p|. The ground state equal time correlations of this latter model are 
known to correspond to the eigenvalue statistics at (any) fixed strength of the perturbation 
for an appropriate choice of the coupling constant, depending on the universality class of the 
perturbation (real, complex or symplectic). If the equivalence between the two systems is 
valid, the full time dependent correlation functions of this model are then known from Ref. 
||T|], and constitute a dramatic progress in our understanding of this many body problem. 
Although the complete spectrum of an integrable many-body problem can often be found, 
the time-dependent correlations involve matrix elements, and the above calculation is one 
of the very few for a non-trivial many body system. 

In this paper, by considering the Hamiltonian of the perturbing potential to be a random 
matrix in an appropriate ensemble 0], we show the equivalence of the evolution of the 
eigenvalues with the dynamics of the Sutherland model. To be more precise, by integrating 
over all perturbations from a given ensemble, we show that the eigenvalues evolve according 
to a Fokker-Planck equation that is equivalent (under a Wick rotation) to the Calogero 
model [^. In the thermodynamic limit, the Calogero model has identical bulk properties as 
the Sutherland model, provided the constants scale properly with N, the size of the matrix 
1^; this scaling is shown to be fulfilled, whereby we establish the conjecture of Ref. [|I|. 

For the evolution of eigenvalues, Simons et al consider a system with a Hamiltonian 
H = Hq + xV, where V is the perturbation, and x is the strength of the perturbation. Units 
are normalized so that the mean level spacing of the eigenvalues of Hq is unity, as is the 
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rms "velocity" of the eigenvalues, defined as {{dei / dx)'^) . The autocorrelation function of 
the energy eigenvalues is defined as 



kix;u;) = (Y,5{e-eiix))5ie-u;-ejix + x)))-l. (1) 

ij 

By explicit calculation for disordered systems, Simons et al provide strong evidence that 
k is universal, with a form that depends only on whether H is in the orthogonal, unitary 
or symplectic ensemble. Numerical simulations are used to argue that this is also true for 
quantum chaotic systems. 

For X = 0, is known to be the same as the ground state equal time correlations in the 
positions of a collection of particles confined to a circle, interacting with the Sutherland 
Hamiltonian. Here we consider the alternative Calogero model in which the particle 
positions are allowed to range over (— oo, oo), with a confining parabolic potential: 

The width of the confining potential is chosen to give a mean ground state inter-particle 
separation in the center of the distribution to be unity. This requirement yields 

= 2N/7r^(3. (3) 

The coupling constant j3 is equal to 1,2 and 4 when H is in the orthogonal, unitary and 
symplectic ensembles respectively. Under the mapping |^ 

= 2it, e = A (4) 

Simons et al argue that k{x; u) is equal to the corresponding time dependent correlation 
function for the Sutherland model. This correlation function can be calculated explicitly for 
(3 = 2, where it agrees with the expression obtained by Simons et al. For (3 = 1 and 4, while 
a complete calculation is not possible for the Sutherland model, the asymptotic forms of k 
for large and small t can be calculated; these agree with the asymptotics of the expressions 
derived for k{x; q) ^ for the orthogonal and symplectic ensembles respectively. 



In this paper we demonstrate that, when is a Gaussian random matrix from the 
appropriate ensemble, the evolution of all quantities such as k, involving only the eigenvalues 
ei{x) at two different values of x, are equal to the corresponding time dependent correlation 
function for the Calogero model. However, correlation functions involving quantities at more 
than two values of x are in general different, so that the full dynamics for the two systems 
are not identical. 

In order to do this, it is more convenient to work with a Fokker-Planck equation that is 
equivalent to the Calogero Hamiltonian. With P{{\i};t) = ipQ{{\i})ip{{Xi};t), where ipo is 
the ground state wave function for the Calogero model, we obtain 



dP{{Xi};t) _^d'P ^ d 



P 



-V— - P +V — 



Xi — Xj 



d 



(5) 



Eq.(|D is then changed to 

x^ = 2r, e = A. (6) 

Eq. is the Fokker-Planck equation corresponding to the Langevin dynamics of a collection 
of classical particles at finite temperature, with logarithmic repulsive pairwise interac- 
tions, the Wigner-Dyson Coulomb gas [|^. Due to the repulsion between the particles, no 
steady state distribution is achieved without the parabolic confining potential. Correlation 
functions that involve the particle positions at two different times ti and T2 > ti can be 
found by considering the evolution from a general initial state for a time T2 — Ti, and then 
averaging over initial states. 

For the quantum system, we modify the parametrization of the disorder strength from 
H{x) = Ho + xV to 

H{x) = Ho cos{nx) + V sm{nx)/n. (7) 

Here V is taken to be of the form V = J2r=o Vr^r, a la Dyson, where Cr are units of the 
appropriate algebra, so that (V^j) = Z]r=o((K)?j) = /^/2(V^i). This parametrization has the 
advantage that when V/Q and Hq are considered to be Gaussian random matrices from the 



same distribution, the distribution for H{x) is stationary as a function of x. For unit mean 
eigenvalue spacing, we obtain 

(Vi) = 2n'N/n'(3, = n'N/n'. (8) 

From the additional normalization condition that the rms velocity of the eigenvalues must 
be unity, it is easy to see from first order perturbation theory that {V^) = 1, so that we 
must choose 



n = ^1x^13 /2N = l/a. (9) 

Since J7 ^ for large A^, the change we have made in parametrizing the disorder strength 
is inconsequential for finite x. 

We first note that, for any perturbation V, the evolution of the eigenvalues e i(^x^ as a 
function of x can be expressed in the form of a set of first order differential equations in the 
continuously changing eigenbasis of H: 

dx " 



/ 12 



dx ti ej 

Eq.(|lOl) can be viewed as a classical Newtonian system, with N degrees of freedom corre- 
sponding to the eigenvalues e^, and N + PN{N — l)/2 degrees corresponding to the diagonal 
and off-diagonal elements of H', obeying appropriate Poisson Bracket relations. This clas- 
sical system can be shown to be integrable [§]; this is because the matrices H and V have 
a very simple dependence on x. The integrability of this system implies that it is essential 
to make assumptions about the nature of the matrix V, either of the form we have made 
above or otherwise. 

It is possible to integrate these equations of motion formally; irreversible equations in 
X arise when Vij are averaged over. In the present work, we find it convenient to follow 



a different strategy to obtain a Fokker-Planclc equation for the eigenvalues. Our method 
follows closely along the lines of a beautiful proof by Dyson |^ for the behaviour of the 
eigenvalues of a matrix subject to random thermal noise; although our problem does not 
have a source of thermal noise, this will be seen to be unimportant for two time correlation 
functions. 

Since V is taken to be a Gaussian random matrix, for any given Hq the matrix H{x) has 
the distribution 

V{H)[dH] = \d{H/ sm{nx))] exp ( - l-Tr\{H{x) - Hocos{nx)}^]—^-—) . (11) 
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We recognize this distribution as the solution to the equation 



sin (ilx) 



dV 



d^V 1 d 



Q,^ H 9^,1^ + ^j^{H^,V) (12) 



with the initial condition H = Hq, where 

aq = 1/Q 

tq = — Oq In [cos(fix)] 

gij = Sij + (1 - Sij)P/2. (13) 

This is the Fokker-Planck equation for a system in which all the matrix elements Hij 
undergo independent Langevin dynamics (starting from inital values (ifo)ij) in a parabolic 
confining well, with an appropriate width and time coordinate. The temperature at which 
the Langevin motion takes place is unity for the diagonal elements of the matrix and (3/2 for 
the offdiagonal elements. Eq.(0) thus implies that, for any given Ho, all moments of the 
different elements of H{x) will be exactly the same at any x as if the elements of H were 
moving independently in parabolic confining wells at a finite temperature. But for such a 
thermal motion, as shown by Dyson 0, it is possible to go to the eigenbasis of H at any 
time, and obtain to second order in perturbation theory for an infinitesmal increase in time 
6tq a simplified form of Eqs.(p!0|): 

6e, = 6H,, + 2 V(e, - e,). (14) 
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The normalization of Eq. ([T3|) then imphes that 



a: 



'Q j^i 



(15) 



The distribution of eigenvalues, PQ{{^i\'-,TQ) at any time tq then satisfies a Fokker-Planck 
equation like Eq.(|^), with Aj replaced by ej, and r and a replaced by tq and aq respectively. 
But comparing Eqs.(|]) to Eqs.® and ([13|) , we see that, in the large N limit, a = aq and 
T = Tq, yielding Eq.®. Thus the distribution of eigenvalues of the quantum system evolves 
in exactly the same way as the distribution of particle positions in a Wigner-Dyson gas. 

Note that we have found that the two distributions evolve in the same way for arbitrary 
initial conditions. Quantities such as k{x;q) can be seen from Eq.(|l]) to involve a sum 
over various moments of e, weighted suitably and then averaged over equilibrium initial 
conditions, requiring only a weaker equivalence. Under a self-averaging assumption, such 
equilibrium averages will be the same as for a generic choice of initial conditions taken from 
the equilibrium distribution. 

Since our results are independent of initial conditions, it is also possible to dispense 
with the parabolic confining well: although there is no longer any steady state, one can 
follow the transient dynamics. This actually corresponds to the original parametrization 
of H = Hq + xV; although conceptually slightly subtle, the algebra is actually simpler. 
Although time translational invariance is now broken, for any fixed finite time the correction 
terms vanish at large A^. 

The result obtained above has been for explicit averaging over a Gaussian random per- 
turbing potential. Apart from the issue of self-averaging, which is relevant for quantum 
chaos, it is necessary to verify that the additional non-Gaussian terms in the distribution of 
V (properly scaled with A^) do not affect the result in the large A^ limit, in order to claim 
universality. We hope to return to this problem in the future. 

While the equivalence between the time coordinate defined in Eq.(|TE|) and Eq.(|^) is 
true only in the large A^ limit, the form given in Eq.(|l^) is true for any N. From the 
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form of Eq.(0), we see that 2Ti/Vt is the Poincare recurrence 'time' interval in x for all the 
eigenvalues to return to their initial values. From Eq.(p^) we see that the interval (0,7r/2) 
in X is stretched out to (0, oo) in tq, so that the system continues to lose memory of its 
initial configuration for all tq, with complete equilibration achieved only in infinite time. 

The ground state correlation function of Eq.(|I|) can be expressed for Gaussian random 
matrices Hq and V in the closed form expression 

^(^i ^) = I ■ wn M / dHoTi \6{e - u - H)] Tr \6{e - Ho) 
I sm [\LX)\ J L J L J 

{ - 2^d^T^'^» + - - 1. (16) 

with u = N + (3N{N — l)/2. All the dependence on time is explicitly present in the cos(nx) 
and sin(f2x) factors. It is sometimes possible to evaluate generalized two-matrix Gaussian 
integrals |]10| , and it would be interesting to apply these techniques to calculating the above. 



It is important to realize that, while Eq.(|TTp implies that all moments of the eigenvalues 
calculated at any x will be equal to the corresponding moments of the particle positions 
of a Wigner-Dyson gas undergoing Brownian motion, this does not mean that the motion 
of the eigenvalues is indeed Brownian. The randomness in the dynamics of the eigenvalues 
comes from the matrix V, which acts like quenched disorder. As a simple illustration of 
the result of the disorder being quenched, we consider the case of = 1. Eq.([TO|) yields 
e(x) — e(0) = V sm{flx)/Q + e(0)[cos(f2a;) — 1]. Averaging over V, for small x we find that 



6(xi) - e{0)Mx2) - e(0)]) = x^X2 + e\0)[Q'xlxl/4]. (17) 

For the Wigner-Dyson gas, on the other hand, for A^ = 1 we have a particle in a parabolic 
well with thermal noise, so that 

[X{n) - A(0)][A(r2) - A(0)]) = 2mm[n,T2] + X'mnT2/a% (18) 



When Ti = T2, with 2r = x^, Eqs.([T7|) and ( p^Sf ) are identical (since a = 1/Q from Eq.(^). 
However, when ti ^ T2, the two equations are different; this difference persists even for large 
A^. 
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It is precisely such multiple time averages that are involved in three point (and higher or- 
der) correlation functions. For instance, the density correlation function that is an extension 
of Eq.(|l|) is 

k{xi,X2;uJi,uj2) = ( ^ ^»(^))'^(^ ^ ^1 - + xi))5{e- UJ2 - ej(x + X2)) ). (19) 

\ ijl I 

This involves various moments of the eigenvalues at two values of x, x\ and X2 (averaged 
over initial conditions). For the Wigner-Dyson gas, the three point correlation functions 
can be expressed in terms of the two point functions, since the first measurement 'rezeroes' 
time. If the two point functions could be calculated for the eigenvalues with arbitrary initial 
conditions {ej(0)}, it would be possible to obtain all higher order correlation functions for 
the Wigner-Dyson gas. Alternatively, even with only a result for a representative choice of 
{ej(0)}, it might be possible to extract equilibrium multiple time correlation functions of 
the Wigner-Dyson gas, if self-averaging is valid. 

In this paper, we have proved the equivalence of the dynamics of the Wigner-Dyson gas 
(or alternatively, the Sutherland model) with the evolution of the eigenvalues of a Hamilto- 
nian under a perturbation drawn from a Gaussian orthogonal ensemble for two time correla- 
tion functions, as argued by Simons, Altshuler et al Multiple time correlation functions 
are not the same in general for the two systems, because of the difference between annealed 
and quenched randomness. It is also possible to prove the equivalence of the two time 
velocity- velocity correlation function, c{u!;x), that measures the correlations in the rate of 
change of the eigenvalues [|I|], by a perturbative method, to all orders in perturbation theory. 
(This method can also be used to obtain the Fokker-Planck equation derived above.) Similar 
considerations should apply to variants of the quantum systems considered here, leading to 
Dyson's circular ensemble j^. The details will be published in a longer paper. 
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